NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


AD-A274  921 


THESIS 


DTIC 


m  ELECTE 


JAN2  5  1994 


A  FINITE  WAKE  THEORY 
FOR  TWO-DIMENSIONAL  ROTARY  WING 
UNSTEADY  AERODYNAMICS 


Thesis  Advisor: 


Mark  A.  Couch 


September,  1993 


E.  Roberts  Wood 


Approved  for  public  release;  distribution  is  unlimited. 


94-02106 

■  _  _ .....  am 


A  A  O 

V/  V 


B4  I  25 


REPORT  DOCUMENTATION  PAGE 

Form  Approved  0MB  No.  0704- 

Public  leportiDg  burden  for  this  coUectioo  of  infotnutioo  is  estinuled  to  average  I  hour  per  response,  including  the  time  for  reviewing  insuucfion,  seardtiog 
existing  data  sources,  gathering  and  truintaining  the  data  needed,  and  coni{Metiog  ai>d  reviewing  the  collectioo  of  information.  Send  comments  regarding  this 
burden  estimate  or  any  other  aspect  of  this  collectioo  of  infatraaiioo,  including  suggestions  for  reducing  this  burden,  to  Wasbingon  headquarters  Services, 
Directorate  for  Information  Operations  and  Reports,  121S  JelTetsoo  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management 
and  Budget,  Paperwork  Reductioo  Project  (0704-0188)  Washington  DC  20S03. 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

23  September  1993  Master’s  Thesis 

4.  TITLE  AND  SUBTITLE  A  RNITE  WAKE  THEORY  FOR  TWO- 
DIMENSIONAL  ROTARY  WING  UNSTEADY  AERODYNAMICS 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S)  Mark  A.  Couch 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey  CA  93943-5000 

8.  PERFORMING 

ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS{ES) 

10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  I 
the  official  policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government,  | 

12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

A 

13.  ABSTRACT  (maximum  200  words) 

The  unsteady  aerodynamic  forces  and  moments  of  an  oscillating  airfoil  for  the  fixed  wing  case  were 
determined  by  Theodorsen  along  with  the  development  of  a  lift  deficiency  function.  Lxiewy  subsequently 
developed  an  analogous  lift  deficiency  function  for  the  rotary  wing  case  in  which  there  are  an  infinite 
number  of  layers  of  shed  vorticity,  or  wakes,  below  the  reference  airfoil.  With  the  advent  of  computer 
panel  codes  that  calculate  the  time  histories  of  the  wakes  generated  by  oscillating  airfoils,  a  theory  is 
developed  for  the  rotary  wing  case  in  which  there  are  a  finite  number  of  layers  of  shed  vorticity  telow 
the  reference  airfoil.  This  theory  includes  a  lift  deficiency  function  that  is  completely  analogous  to 
Loewy  and  Theodorsen. 

It  has  long  been  recognized  that  an  airfoil  oscillating  in  pure  plunge  produces  a  propulsive  force 
("Katzmayr  effect").  Garrick  used  Theodorsen's  work  to  develop  equations  for  the  propulsive  force  that 
include  the  lift  deficiency  function  as  a  parameter.  When  either  Loewy's  lift  deficiency  function  or  the 
finite  wake  lift  deficiency  function  is  used,  the  effect  of  the  propulsive  force  is  greatly  enhanced  with  the 
proper  phase  relationship  of  the  wakes.  The  finite  wake  theory  along  with  Garrick's  work  is  used  to 
describe  the  performance  characteristics  of  Higher  Harmonic  Control.  Specifically  for  the  OH-6A, 
coupled  pitch-plunge  motion  results  in  a  propulsive  force  that  significantly  reduces  the  rotor  drag  force. 


14  SUBJECT  TERMS 

Helicopter,  Rotary  Wing,  Unsteady  Aerodynamics,  Wake,  Propulsive  Force, 

Higher  Harmonic  Control,  Hummingbird 

15.  NUMBER  OF 

PAGES  69 

16.  PRICE  CODE 

17.  SECURITY  CLASSin- 
CATION  OF  REPORT 
Unclassified 

18.  SECURITY  CLASSIFI¬ 
CATION  OF  THIS  PAGE 
Unclassified 

19.  SECURITY  CLASSIFICA¬ 
TION  OF  ABSTRACT 
Unclassified 

20.  UMITATION  OF 
ABSTRACT 

UL 

NSW  7540-0 1-280-3500  Standard  h'orin  298  tKev.  2-89) 


Prescribed  by  ANSI  Std.  239-18 


1 


lyno  QUALITY  INSPECTED  6 


Approved  for  public  release;  distribution  is  unlimited. 

A  Finite  Wake  Theory 
For  Two-Dimensional  Rotary  Wing 
Unsteady  Aerodynamics 

by 

Mark  A.  Couch 

Lieutenant,  United  States  Navy 
B.S.,  Purdue  University,  1984 

Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  AERONAUTICAL  ENGINEERING 

from  the 


Author: 


NAVAL  POSTGRADUATE  SCHOOL 


September  1993 


Mark  A.  Couch 


Approved  by: 


Max  F.  Platzer,  Second  Reader 

Daniel  J.  Chairman 

Department  of  Aeronautics  and  Astronautics 


ii 


ABSTRACT 


The  unsteady  aerodynamic  forces  and  moments  of  an  oscillating  airfoil  for  the 
fixed  wing  case  were  determined  by  Theodorsen  along  with  the  development  of  a  lift 
deficiency  function.  Loewy  subsequently  developed  an  analogous  lift  deficiency 
function  for  the  rotary  wing  case  in  wliich  there  are  an  infinite  number  of  layers  ot 
shed  vorticity,  or  wakes,  below  the  reference  airfoil.  With  the  advent  of  computer 
panel  codes  that  calculate  the  time  histories  of  the  wakes  generated  by  oscillating 
airfoils,  a  theory  is  developed  for  the  rotary  wing  case  in  which  there  are  a  finite 
number  of  layers  of  shed  vorticity  below  the  reference  airfoil.  This  theory  includes  a 
lift  deficiency  function  that  is  completely  analogous  to  Loewy  and  Theodorsen. 

It  has  long  been  recognized  that  an  airfoil  oscillating  in  pure  plimge  produces  a 
propulsive  force  ("Katzmayr  effect").  Garrick  used  Theodorsen's  work  to  develop 
equations  for  the  propulsive  force  that  include  the  lift  deficiency  function  as  a 
parameter.  When  either  Loewy's  lift  deficiency  function  or  the  finite  wake  lift 
deficiency  function  is  used,  the  effect  of  the  propulsive  force  is  greatly  enhanced  with 
the  proper  phase  relationship  of  the  wakes.  The  finite  wake  theory  along  with 
Garrick's  work  is  used  to  describe  the  performance  characteristics  of  Higher  Harmonic 
Control.  Specifically  for  the  OH-6  A,  coupled  pitch-plunge  motion  results  in  a 
propulsive  force  that  significantly  reduces  the  rotor  drag  force. 
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TABI  i:  t  F  SYMBOLS 


a  non-dimensional  elastic  axis  1<  .tic  i  measured  from  midchord 
A  mathematical  expression  dencti  g  g  oup  of  terms  described  by  equation  (61) 

An  mathematical  expression  dencti  g  g  oup  of  terms  described  by  equation  (77) 

b  semichord  length 

B  mathematical  expression  dencti  .g  g  oup  of  terms  described  by  equation  (62) 

Bn  mathematical  expression  denciing  g  oup  of  terms  described  by  equation  (78) 

C(k)  Theodorsen's  lift  deficiency  fut.ctio  .,  C(k)  =  F  +  iG 
C'(k)  Loewy's  lift  deficiency  function.,  C'(k)  =  F'  +  iG’ 

C*(k)  finite  wake  lift  deficiency  function,  C'(k)  =  F"  +  iG' 

Cpj  propulsive  force  coefficient 
F  real  part  of  Theodorsen's  lift  deficiency  function 
F'  real  part  of  Loewy's  lift  deficiency  function 
F*  real  part  of  finite  wake  lift  deficiency  function 
G  imaginary  part  of  Theodorsen's  lift  deficiency  function 

G'  imaginary  part  of  Loewy's  lift  deficiency  function 
G*  imaginary  part  of  finite  wake  lift  deficiency  function 
h  non-dimensional  distance  betv  een  layers  of  shed  vorticity  (wake  spacing) 
h  - 

"  bQCl 

b  simple  harmonic  motion  of  venlca!  deflection,  h  = 
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ho  amplitude  of  deflection  in  simple  harmonic  motion  of  (plunge) 

ho  non-dimensional  amplitude  of  deflection,  h.  =  — 

b 

h  first  derivative  with  respect  to  time  of  vertical  deflection  motion  (plunge) 
h  second  derivative  with  respect  to  time  of  vertical  deflection  motion  (plunge) 

Ho  Hankel  function  (complex  Bessel  function)  of  second  kind  of  order  n, 

i  complex  notation  ,  i  = 

Jo  real  part  of  complex  Bessel  function  of  order  n 
k  reduced  frequency,  k  = 

V 

m  ratio  of  oscillation  frequency  to  rotational  frequency,  tn  = 

«  it 

Mo  aerodynamic  moment  about  the  elastic  axis  (positive  clockwise) 

n  revolution  number 

N  number  of  rotor  revolutions  (reduces  to  the  number  of  wakes  for  single-blade 
rotor) 

Ap  pressure  difference  across  airfoil 

P,  propulsive  force 

q  blade  number 

Q  total  number  of  rotor  blades 

r  blade  section  radius 

S  x-component  of  induced  downwash  velocity 

t  time 

U  vortex  distribution  per  unit  length 
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V  freestream  velcxiity 

V  forward  airspeed 

w  downwash  velocity 

W  Loewy's  wake  weighting  function 
Wn  finite  wake  weighting  function 

X  non-dimensional  distance  from  midchord 

Xo  non-dimensional  distance  from  midchord 

Yn  imaginary  part  of  complex  Bessel  function 
a  simple  harmonic  motion  of  pitch  angle,  a  = 
tto  amplitude  of  rotation  in  simple  harmonic  motion  of  aoe'“‘ 
d  first  derivative  with  respect  to  time  of  rotation  motion  (pitch) 

&  second  derivative  with  respect  to  time  of  rotation  motion  (pitch) 
d  real  part  of  Loewy's  wake  weighting  function 
On  real  part  of  finite  wake  weighting  function 
^  imaginary  part  of  Loewy's  wake  weighting  function 
imaginary  part  of  finite  wake  weighting  function 
Y,  vorticity  generated  by  reference  airfoil 
Yn,  vorticity  generated  by  9th  blade  in  mh  revolution 
r  total  circulation  around  airfoil 

p  advance  ratio,  p  = 

fir 

5  non-dimensional  distance  from  midchord 

p  density  of  air 
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I.  INTRODUCTION 


A.  GENERAL 

The  present  knowledge  of  rotary  wing  unsteady  aerodynamics  is  limited.  The 
fundamental  closed  form  solutions  by  Theodorsen  [Ref.  1]  and  Loewy  [Ref.  2]  provide 
the  basis  for  the  theoretical  work  in  this  area.  While  Loewy's  work  on  wake-induced 
flutter  helps  explain  the  phenomenon,  it  also  points  out  the  difficulty  to  be  overcome 
which  is  the  closed  form  theory  shows  rapid  changes  in  the  lift  deficiency  function 
with  changes  in  reduced  frequency  (k),  wake  spacing  {h),  and  frequency  ratio  (w). 

In  the  past,  emphasis  in  the  study  of  unsteady  aerodynamics  has  tended  to  focus  on 
flutter  instability  and  the  effects  of  unsteady  aerodynamics  on  lift  generation  and 
torsional  loads,  and  not  on  performance.  In  this  thesis,  the  emphasis  will  be  on 
performance  and  the  effect  of  unsteady  aerodynamics  on  the  drag  of  the  airfoil.  The 
classic  work  on  this  subject  was  performed  by  Garrick  [Ref.  3].  In  this  work,  Garrick 
found  hat  the  drag  force  could  be  reversed,  acting  to  propel  the  airfoil  forward.  This 
drag  reversal  is  generally  referred  to  as  either  negative  drag  or  propulsive  force. 

The  primary  objective  of  this  thesis  is  twofold.  First,  the  effect  of  rotary  wing 
unsteady  aerodynamics  on  the  drag  of  an  airfoil  will  be  analyzed.  It  will  be  shown 
that  this  effect  can  be  either  substantially  greater  than  that  shown  by  Garrick,  or 
substantially  less  depending  on  the  frequency  ratio  {m).  This  change  in  the  effect  is 
due  to  layers  of  shed  vorticity  (or  wakes)  below  the  rotor  blades.  Second,  with  the 
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recent  advent  of  computer  panel  code  c:  >able  of  calculating  the  time  histories  of  the 
wakes  generated  by  oscillating  airfoiU.  a  dosed  form  theory  is  developed  for  the 
rotary  wing  aerodynamics  in  which  t.  re  at  a  finite  number  of  wakes  below  the  rotor 
blades.  It  will  be  shown  that  this  th.  ry  is  completely  analogous  to  Theodorsen  and 
Loewy's  work,  and  applicable  to  Gar  .k',  work. 

B.  SCOPE 

Chapter  D  is  a  short  review  of  th;  wc  k  of  Theodorsen  [Ref.  1]  and  Garrick  [Ref. 
3].  The  material  is  presented  in  a  n:  me  to  allow  one  to  understand  the  key 
assumptions  made  in  each  theory  v  iil.out  being  overcome  by  the  details.  Chapter  ni 
contains  an  analysis  of  Loewy's  work  [Rt  f.  2]  and  applies  it  to  Garrick.  The  closed 
form  finite  wake  theory  is  also  presented  here,  as  well  as  the  special  case  of  a  single 
wake.  Chapter  IV  is  a  limited  discussion  on  applications  of  finite  wake  rotary  wing 
unsteady  aerodynamic  theory.  The  fust  section  discusses  logical  extensions  of  this 
theory  to  include  forward  airspeed  [Re‘.  4]  and  compressibility  effects  [Ref.  5  and  6]. 
A  comparison  of  this  theory  to  the  Na  va!  Postgraduate  School  (NPS)  Unsteady  Panel 
Code  [Ref.  7  and  8]  for  the  case  of  pure  plunge  is  also  included  in  the  first  section. 
The  second  section  presents  an  explaiiation  for  the  performance  benefits  due  to  Higher 
Harmonic  Control  (HHC).  Appendix  A  is  another  application  of  this  theory  describing 
the  efficient  nature  of  flight  of  the  huinrungbird.  The  discussions  in  Chapter  IV  and 
Appendix  A  are  not  designed  to  provide  all  the  answers  to  the  questions  arising  from 
these  applications.  They  are  designed  tc  show  that  positive  performance  effects  can  be 
found  in  the  application  of  rotary  wing  i  nsteady  aerodynamics. 
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n.  BACKGROUND 


A.  THEODORSEN'S  LIFT  DEFICIENCY  FUNCTION 

Tlie  aerodynamic  forces  and  moments  on  an  oscillating  airfoil  for  the  fixed-wing 
case  as  determined  by  Theodorsen  [Ref.  1]  are  based  on  potential  flow  theory  and  the 
Kutta  condition.  Potential  flow  theory  reduces  to  Laplace's  equation  for  the  case  of 
incompressible  subsonic  flow. 


=  0 


(1) 


dx^  dy^ 

Since  potential  flow  theory  follows  the  principle  of  superposition,  the  potentials  can  be 
treated  as  two  types;  non-circulatory  and  circulatory.  The  non-circulatory  terms  are 
primarily  from  the  effects  of  the  freestream  with  the  airfoil.  The  circulatory  terms 
result  from  the  vorticity  generated  in  the  wake  of  the  airfoil.  It  is  in  the  circulatory 
terms  that  the  lift  deficiency  function  is  defined. 

The  aerodynamic  forces  and  moments  on  a  thin  airfoil,  depicted  in  Figure  1,  are 
obtained  by  integrating  the  pressure  difference  across  the  airfoil  determined  from  the 
generalized  Bernoulli  equation  as  follows: 


P  = 


p{x)dx 


(2) 
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T* 

-/ 

trailing 

edge 

All  distances  are  non-dimensionalized  by  the  semi-chord  (b) 

Figure  1.  Parameters  of  the  thin  airfoil. 


and 

j  p{x){x-  a)dx  ,  (3) 

where  the  generalized  Bernoulli  equation  is  written  as 

.pw  =  -2p(^i‘E  .  .  (4) 

'  b  cx  dt ' 

In  the  manner  suggested  by  Scanlan  and  Rosenbaum  [Ref.  9],  the  potential  flow 
for  the  non-circulatory  terms  can  be  written  as 
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CP.  =  ^log 


(5) 


TTie  pressure  difference  across  the  airfoil  for  the  non-circulatory  terms  is 

,ryr?J-.  •'-c  (6, 

- f  ‘uizJL:^£ISUa)ai  , 

^  Ji  4-;r5  +  \/T^/rF 

where  the  downwash*  w(^)  =  h  +  b(i-a)d  +  va,  and  simple  harmonic  motion  is 
assumed,  or  tv  =  Woe'"'.  Integrating  the  pressure  difference  in  equation  (6)  to  obtain 
the  non-circulatory  forces  and  moments  yields 

/*.  =  -TTfih^ivd +£- baa) 


and 

(A^).  =  -iTpA^|(^-a)KM  +  A^(-|+a^)d: -flAirj  .  (8) 

The  potentials  of  the  circulatory  terms  are  determined  by  integrating  the  velocity 
potential  resulting  from  a  single  vortex  element  with  a  vortex  distribution  per  unit 
length,  U(xJ,  that  is  assumed  to  vary  sinusoidally.  This  yields  a  circulation  potential 
of 


‘Only  two  degrees  of  freedom  (pitch  and  plunge)  are  presented  here  since  most  rotary 
wing  aircraft  do  not  have  ailerons  incorporated  on  the  rotor  blades. 
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-1  [ 

— - - tan  - - - 

2ir  I- x-x^ 


\dx„  . 


The  Kutta  condition,  which  requires  the  velocity  at  the  trailing  edge  to  be  finite,  can 


be  written  as 


'  OX  'edge  V  dx  'edge 


This  results  in 


1-?  2j, 


Now  defining  a  lift  deficiency  function  as 


ak)  = 


r-^e-“ 

■'<  i/^TT 


the  Kutta  condition  becomes 


-2wa)dL 


Him)  = 


Using  equation  (4)  the  pressure  difference  across  the  airfoil  for  the  circulatory  terms  is 
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t^p{x)Y  = 


'TrV'l  -X 


ill 
N  1-5 


(14) 


Integiatmg  the  pressure  difference  in  equation  (14)  to  obtain  the  circulatory  forces  and 
moments  yields 


/>r  =  -2npvbak)[va-Ji-b{~-a)a] 


(15) 


and 


(/K.)r  =  2irpvi’(a*-i)ai)[vn*irf'(|-aW  ■ 


(16) 


The  total  pressure  difference  is  obtained  by  adding  the  pressure  differences  for  the 
non-circulatory  and  circulatory  terms.  After  some  algebraic  manipulation  the  total 
unsteady  pressure  distribution  can  be  written  as 


i^p(.x) 

pv 


yx 


1-5  1 


2 


N  i-5V^-5 

1  'Xl  -v/l  -X^^\  -5^ 


(17) 


w{l)dl 


Likewise  the  total  forces  and  moments  are  obtained  by  adding  their  two  parts  together 
yielding 


P  =  -TTpb'^(va*h'baa)-27rpvbak-)\va^h^b[Y^)a\ 


(18) 
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and 


+  27r p  . .  +  -  -)  (Xk) I  m  -^  ^  +  ^ (^  -  a)dj  . 

The  analytical  solution  to  the  lif;  efi  iency  function  involves  the  use  of  complex 
Bessel  functions  (or  Hankel  functio:.^  ,  a?  d  can  be  written  as 


wwW 

ak)  55  . 


where  -  Jn-  *^i  ts  the  Hankel  funct!  jn  of  the  second  kind  of  order  n.  The  lift 
deficiency  function  can  be  separated  into  its  real  and  imaginary  parts  as 


CXk)  -  F{k)*iG{k)  , 


where 


m  = 


/,(7|.  Ko)«  r|(yi-/„) 
(7, .  ♦  ( r,  -7„)» 


Kr„«7,7„ 

CKk)  - -!-2 — — 7  , 

(7, .  r„)>  *  ( r,  -7o)» 


and  all  Bessel  functions  are  evaluated  at  the  reduced  frequency  {k).  A  typical  plot  of 
F(k)  and  G(k)  as  given  by  Theodorsen  [l!ef.  1],  is  shown  in  Figure  2.  Figure  3  is  a 
semi-logarithmic  plot  of  Theodorsen's  lif  deficiency  fimction,  and  the  limits  of  F(k) 
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F  or  -G 


Figure  2.  Conventional  plot  of  Theodorsen's  lift  deficiency  function. 


Figure  3.  Semi-logarithmic  plot  of  Theodorsen's  lift  deficiency  function. 
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and  G(k)  as  k  approaches  zero  and  infinity  can  be  more  easily  seen.  It  is  an 

interesting  note  that  when  =  0,  t  =  -G  *  0.188773655.,.. 

ok 

B.  GARRICK'S  PROPULSIVE  FORCE  EQUATION 

One  of  the  earliest  investigations  on  the  effect  of  an  oscillating  airfoil  on  the 
resultant  drag  was  conducted  by  KatzmajT  [Ref.  10]  in  two  series  of  experiments.  In 
his  second  series,  Katzmayr  measured  a  negative  drag  or  propulsive  force  for  certain 
velocities  and  frequencies.  This  discovery,  later  to  be  known  as  the  "Katzmayr 
effect",  helped  answer  questions  on  how  birds  propel  themselves  through  the  air. 
Richardson  [Ref.  11]  used  the  "Katzmayr  effect"  to  explain  the  locomotion  of  fish 
through  the  water,  Garrick  [Ref.  3]  applied  the  equations  of  motion  determined  by 
Theodorsen  [Ref,  1],  equations  (18)  and  (19),  to  a  method  developed  by  von  Karman 
and  Burgers  [Ref.  12]  to  obtain  a  closed  form  analytical  solution  for  the  propulsive 
force  generated  by  an  oscillating  airfoil.^ 

Garrick  determines  the  average  horizontal  force  by  two  methods:  energy  and 
force.  He  then  compares  the  results  of  these  methods,  and  shows  they  are  identical. 
For  simplicity  only  the  force  method  is  presented  here.  First  he  assumes  simple 
harmonic  motion  for  each  of  the  degrees  of  freedom,  resulting  in 

a  =  (24) 


^Though  Garrick  developed  his  equations  for  three  degrees  of  freedom,  only  two 
degrees  of  freedom  (pitch  and  plunge)  will  be  presented  here. 
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and 


h  =  .  (25) 

where  (pg  and  (pj  are  phase  angles  between  the  initiation  of  the  input  and  an  arbitrary 
reference  point. 

The  total  horizontal  force  on  the  airfoil  is  given  by 

Pjf  =  TTpS^^OtP  (26) 


where  P,  is  the  propulsive  force  on  the  airfoil  (positive  to  the  left),  P  is  the  resultant 
force  on  the  airfoil  (positive  down),  and  S  is  the  component  of  induced  downwash  in 
the  x-direction.  The  time  averaged  propulsive  force  is  found  by 

P  .  dt  .  (27) 

2irio 

The  x-component  of  the  induced  downwash  is 


=  — 12 C(ir)|  va+i*  b[^  -  a)dj  -  Z>d|  . 


(28) 


In  finding  the  solution  to  the  time  averaged  propulsive  force  equation,  it  is 
arbitrary  whether  one  chooses  to  employ  the  real  or  imaginary  parts;  the  results  will  be 
the  same.  Ganick  [Ref.  3]  chose  the  imaginary  parts,  and  determined  the  time 
averaged  propulsive  force  to  be 
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(29) 


Pj^  =  irpAw^jA^ 

■"  ^  ^  "  "f ) 

^  ([{  -  4^^  ■^^^)*\{\--J~^))  C0S((P2  -  ‘Po)]}  • 

It  is  interesting  to  note  that  in  the  special  cases  of  only  one  degree  of  freedom,  the 
effects  of  the  arbitrary  phase  angles  between  initiation  and  the  reference  point 
disappear.  The  propulsive  force  equations  become  for  pure  plunge  (a^  =  0) 


F,  =  Trpb(o^/il(F^^G^)  , 


and  for  pme  pitch  (h^  *  0)^ 


'2  '2  >  ki 


If  it  is  assumed  that  pitch  and  plunge  are  in  phase  with  each  other  (they  reach  their 
maximum  deflections  at  the  same  time),  the  coupled  propulsive  force  equation 
becomes 


’The  result  presented  here  differs  from  Garrick's  equation  (34)  in  reference  3.  When 
equation  (29)  in  reference  3  is  simplified  using  the  definitions  of  and  bj,  the  result  is 
equation  (31)  above.  This  result  has  been  coi^irmed  by  Garrick's  work  in  reference  10. 
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=  irpA(i>* 


Since  propulsive  force  is  analogous  to  negative  drag,  it  is  convenient  to  define  a 
propulsive  force  coefficient  as 


P  P 

C  -  ^  =  X 

^pv'^O.b) 


The  propulsive  force  coefficients  for  the  special  cases  become  for  pure  plunge 
(tto  -  0) 


Cp^  =  7rk^hl[F'^*G'^)  ,  (34) 

where  lio  is  the  amplitude  of  deflection  non-dimensionalized  by  the  semi-chord,  for 
pure  pitch  (h^  =0) 

Cp.  ‘  "If}  • 


and  for  coupled  pitch  and  plunge 


Cp,  =  nk^[P,{P^-G^yal[(F^^G^){±*{^-afy^(^^ 
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ID  /  NALYSIS 


A.  LOEWY'S  LIFT  DEHCIENC'  FL  ECTION 

The  aerodynamic  forces  and  moi  nt^  on  an  oscillating  airfoil  for  the  rotary-wing 
case  are  more  complex  than  the  fixe .\  i’  ;  counterpart.  Loewy  [Ref.  2]  determines 
these  forces  and  moments  for  the  ht  \  r  c  se  where  k  =  ,  but  uses  a  slightly 

different  approach  to  the  problem  by  '  jh  ig  the  integral  downwash  equation  which 
leads  to  an  equation  for  the  pressure  distr  jution  which  is  in  the  same  form  as 
Theodorsen.  Loewy  then  defmes  a  n  •  dif  ed  lift  deficiency  function,  and  states  that 
the  integration  of  the  pressure  distribii  ior  across  the  airfoil  would  be  identical  to 
Theodorsen.  This  makes  the  expressions  for  forces  and  moments  identical  except  for  a 
modified  lift  deficiency  ftmetion. 

The  first  step  that  Loewy  perfomi  •  is  to  set  up  a  system  to  account  for  the  wakes 
generated  by  the  previous  blades  in  tho  same  revolution  and  all  blades  in  previous 
revolutions  as  shown  in  Figure  4.  He  u'-es  two  indices  to  account  for  the  vorticity 
shed  by  a  given  wake:  n  indicates  the  revolution  and  q  indicates  the  blade  whose 
wake  it  is.  The  induced  velocity  or  dowii  wash  resulting  from  an  element  of  vorticity 
is  obtained  from  the  Biot-Savart  Theou  iri, 


divix)  = 


2‘7r[U  -  0^  + 
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Figure  4.  Aerodynamic  model  for  a  multi-blade  rotor  system. 


where  is  the  vorticity,  Q  is  the  number  of  rotor  blades,  and  h  is  the  non- 
dimensional  wake  spacing.  Writing  the  integrals  involving  the  bound  vorticity  and  the 
vorticity  in  the  wake  of  the  reference  airfoil  separately  from  the  rows  of  vorticity 
below  the  plane  of  the  rotor  yield 

^(;r)  = 

r  ] 

h\  J-oc  {x-if  *  n'^Q'^h'^l 

In  essence  the  problem  has  been  broken  down  in  the  same  maimer  as  Theodorsen. 
The  first  integral  represents  the  effects  of  the  freestream  on  the  airfoil  (non-circulatory 
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terms).  The  second  integral  represents  the  velocity  created  by  the  vorticity  generated 
by  the  reference  wake  (circulatory  term).  The  third  and  fourth  integrals  represent  the 
velocity  created  by  the  vorticity  generated  by  previous  blades  or  in  previous 
revolutions  (circulatory  terms).  The  main  difference  between  Loewy  and  Theodorsen 
is  the  terms  which  account  for  the  vorticity  generated  by  previous  blades  or  in 
previous  n  volutions. 

Loewy  [Ref.  2]  shows  that  the  vorticity  shed  by  the  qth  blade  in  the  nth  revolution 
is 

Y  =  y^r  (39) 


where  F  is  the  total  circulation  around  the  airfoil,  tj;,  is  the  phase  angle  by  which  the 
motion  of  the  ^h  blade  leads  that  of  the  reference  blade,  and  m  is  the  ratio  of 
oscillatory  frequency  to  rotational  frequency.  Substituting  the  vorticity  expression 
from  equation  (39)  into  the  integral  do^^’nwash  equation  (38)  yields 


ffix)  =  - 


1  1  f '  ,, 


^If 


-  j: 


^-iO.iraqlQ-’h  ) 


a-}  J. 


-nnma  f  g  (XT- $)(/$ 

J-*  (x- 0^  ^ 


The  last  two  integrals  in  equation  (40)  have  the  form 


-«c  {x-lf^A'^ 


(41) 
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Substituting  equation  (41)  into  equation  (40)  and  noting  that  the  sununations  over  n 
are  convergent  geometric  series  give 


2rrlJ_j  (jf-0  J,  (x-0 


(42) 


where 


<?-i 


1  ^  E 

W{kji^)  =  - aii - 

gkbQgilTrm  _  j 


(43) 


The  function  W(k,h,ni}  may  be  thought  of  as  a  weighting  function  for  the  vorticity 
shed  by  previous  blades  or  in  previous  revolutions. 

The  form  of  the  downwash  equation  in  equation  (42)  can  be  solved  by  applying 
Sohngen's  inversion  formula  [Ref.  14],  which  shows  that  the  solution  to  an  equation  in 
the  form 


g^x)  = 


mdi 

X-l 


(44) 


IS 


'  -u 


r~i —  c  * 
1  -X 


Nl-?  .*-5 


(45) 


Satisfying  the  condition /(i)  =  finite  is  the  same  as  employing  the  Kutta  condition. 
The  bound  vorticity  becomes 
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Y,U)  =  - 


1  li£  f 

7r\  1  IJ.j'Nl  1-5  x-i 


‘  n^/  1  \re-“^dr)  .. 


Evaluating  the  circulation  over  the  entire  airfoil  yields 


r  =  f  \/x)dx  =  _ LillZ^ _  ,  (47) 

J-i  '  +  (/,  ^  //o)  ’ 

where  the  Hankel  and  Bessel  functions  are  evaluated  at  reduced  frequency  (i). 

Since  the  airfoil  can  be  thought  of  as  a  vortex  sheet,  the  generalized  Bernoulli 
equation  (4)  becomes 


2/". 

J-jN  1  ~  C 


Aj^(x)  =  -pjyyM/) 

Pij-i 


If  simple  harmonic  motion  is  assumed. 


and  equation  (48)  becomes 


Ap(x)  =  -p  vY,(>*r,r)  + 


Substituting  the  bound  vorticity  equation  (46)  and  the  airfoil  circulation  equation  (47) 
into  equation  (50)  yield 


18 


^pjx) 

pftr 


r 

+  2[/,  +  //o]  Wikji^)  J -1 

■ilM 


1  -X 


_  ijLl 

I+a-N  1-5 


-_£ 

jr 


1  +  5/  1 


(51) 


2 


N  1-5W-5 

1  -jr5  -/r^v/i~^ 

1  -jr5  +^1  -Ar^V^l  -5^ 


logj 


w{l)dl 


The  pressure  distribution  in  equation  (51)  has  the  same  form  as  Theodorsen's  in 
equation  (17)  if  the  factor  multiplying  the  first  integral  is  written  as 


^\\  -  C'{kJijn)\  . 

tj-l 


(52) 


where  C'(k,h,m)  is  Loewy's  lift  deficiency  ftmction.  Solving  for  C'(k,h,m)  yields 


C'i.kJijm') 


^  2/,  WjkJijn) 

/f®  -  ^  2[/,  -  /7o]  W{kji^)  ’ 


(53) 


where  the  Hankel  and  Bessel  functions  are  evaluated  at  reduced  frequency  (it).  Since 
C’(k,h,m)  is  not  a  function  of  location  along  the  airfoil,  the  integration  of  the  pressure 
distribution  in  equation  (51)  across  the  airfoil  will  yield  equations  of  motion  identical 
to  Theodorsen's  in  equations  (18)  and  (19)  except  C(k)  will  be  replaced  by  C'(k,h,m). 

It  can  be  seen  that  as  W(k,h,m)  approaches  zero,  C'(k,h,m)  =  C(k),  and  this  condition 
corresponds  to  an  infinite  wake  spacing  (A-*®). 

Since  the  wake  weighting  function  is  periodic  in  nature,  Loewy  [Ref.  2]  shows  that 
the  wake  weighting  function  for  a  multi-blade  rotor  can  be  expressed  by  a  single-blade 
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20 


(59) 


0>  =  (yi-^A^)A*MU2a)B 
A^^B^ 

and 


(60) 


A  =  J,(U2&)*  y,- 2  Jj  («) 

B  =  -  r,  +  2/,/5  +  J^il  +  2a)  .  (62) 

Plots  of  Loewy's  lift  deficiency  function  are  shown  in  Figures  5  through  8.  It  is  easily 
seen  that  there  are  rapid  changes  in  the  lift  deficiency  function  as  reduced  frequency 
(Jl:),  wake  spacing  (h),  and  frequency  (m)  vary. 

Since  Loewy's  lift  deficiency  function  is  completely  analogous  to  Theodorsen's  lift 
deficiency  function,  the  propulsive  force  determined  by  Garrick  [Ref.  3]  and  the 
propulsive  force  coefficient  can  be  calculated  for  the  rotary-wing  case  by  substituting 
F'  and  G'  for  F  and  G.  For  pure  plunge  (Oo  *  0) 

,  (63) 


for  pure  pitch  (h,,  =  0) 


■/ 


(64) 
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Figure  6.  Loewy's  lift  deficiency  function  as  a  function  of  wake  spacing  for 
m  -  0.25. 
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Figure  7.  Loewy's  lift  deficiency  function  as  a  funciton  of  wake  spacing  for 
m  -  0.5. 


Figure  8.  Loewy's  lift  dieficiency  function  as  a  function  of  wake  spacing  for 
m  -  0.75. 
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and  for  coupled  pitch  and  plunge 


All-,  iF){l-s.±Ul.s)i^ 
2'2  h  P’  ^2  >  k 


(65) 


Plots  of  propulsive  force  coefficient  ;  a  i  inction  of  frequency  ratio  for  k  =  0.1234  are 
shown  in  Figure  9  through  11.  It  c;  i  bo  :  een  that  Cp^  is  reasonably  constant  for  0.15 
<  m  <  0.8  and  greater  that  the  "Katznoivr  effect."  (In  this  case  h  -  100  is  sufficient  to 
represent  For  m  <  0.15  and  m  >  0.;^,  Cp^  begins  to  decrease  to  the  point  where 

it  eventually  becomes  less  than  the  "Ki-tzinayr  effect." 

B.  nNITE  WAKE  LIFT  DEFICIENCY  FUNCTION 
1.  General 

Loewy's  theory  can  be  extended  to  the  case  of  a  finite  number  of  wakes 
below  the  rotor  blade  by  modifying  the  uake  weighting  function.  This  modified  wake 
weighting  function  is  completely  analc^oes  with  Loewy's  which  makes  the  definition 
of  the  lift  deficiency  function  the  same,  b  it  with  a  and  ^  modified  by  the  new  wake 
weighting  function.  Since  the  lift  deficie;  cy  function  is  of  the  same  form,  the 
aerodynamic  forces  and  moments  will  be  ‘he  same,  but  with  a  modified  lift  deficiency 
function.  Loewy  has  already  shown  liiat  he  rotary-wing  problem  reduces  to  the  case 
of  a  single  rotor.  'The  modified  wake  we  ^hting  function  will  also  be  developed  for 
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Figure  10.  Propulsive  force  coefficient  in  pure  pitch  as  a  function  of  wake  pacing 
with  an  infinite  number  of  wakes  (Loewy's  lift  deficiency  function). 
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Figure  11.  Propulsive  force  coefficient  in  coupled  pitch-plunge  as  a  function  of  wake 
spacing  with  an  infinite  number  of  wakes  (Loewy's  lift  deficiency  function). 


the  case  of  the  multi-blade  rotor,  and  then  reduced  to  the  case  of  the  single-blade  rotor 
with  modified  wake  spacing  and  frequency  ratio. 

The  development  of  a  finite  wake  lift  deficiency  function  is  the  same  as 
Loewy's  described  earlier.  Changing  the  integral  downwash  equation  (40)  to  a  finite 
number  (AO  of  revolutions  yields 


wix)  =  - 


1  yprf^JL!^ 

2irlJ.,  j, 

-  /A-rV  g-ilnao  ^  C (.T  -  Q 


(66) 
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where  the  summations  from  n  equal  zero  to  infinity  have  been  replaced  by  finite 
summations  from  n  equal  zero  to  N.  Using  equation  (41)  to  solve  the  last  two 
integrals  in  equation  (66)  results  in 


di 


,0-^  .  ^ 

g-\ 

N  : 

fl-1  ■ 


^-5 
( 


(67) 


Since  this  is  in  the  same  form  as  equation  (42),  the  terms  enclosed  by  the  parentheses 
can  be  defined  as  the  finite  wake  weighting  function,  or 

<p-i  N  N 

^=1  iJ=0  11=1 

It  is  easily  seen  that  this  modified  wake  weighting  function  is  periodic  once  the 
relationship  is  known,  and  a  multi-blade  rotor  can  be  reduced  to  a  single-blade 
rotor  with  modified  wake  spacing  and  frequency  ratio.  For  the  case  of  the  single¬ 
blade  rotor  (Q  =  1), 


N 

^-ilwma ^-Dkb 
11=  I 


(69) 


It  can  be  shown  that  as  N-><»,  the  finite  wake  weighting  function  in  equation  (69) 
becomes  a  convergent  geometric  series  equivalent  to  Loewy's  wake  weighting  function 
in  equation  (54)  provided  m  *  Q  and  h  *  Q,  k  singularity  in  the  solution  exists  when 
m  ■  0  and  A  ■  0,  and  the  case  of  h  *  0  has  no  real  physical  significance. 
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The  solution  to  the  remainder  of  the  finite  wake  problem  is  identical  to 
Loewy,  but  using  the  modified  finite  wake  lift  deficiency  function.  To  keep  the  finite 
wake  lift  deficiency  function  separate  fi-om  Loewy's  lift  deficiency  function,  the  former 


is  defined  as  follows: 


C’ikJijn) 


*  IH^  »  2[y,  *  //„]  ' 


(70) 


Separating  the  finite  wake  weighting  function  in  equation  (69)  into  real  and  imaginary 


parts, 


'N 


(71) 


yields 

N 

a^\ 

and 


N 

(-e'®^*sin  27r/n/7)  . 

11=  1 

The  finite  wake  lift  deficiency  function  becomes 


C*  =  F'^iG'  , 


(72) 


(73) 


(74) 


where 


30 


(75) 


and 

/I;,. /,(1.2d^«  (77) 

B„  •  -r,>2J,0^*J^{\*2a,)  .  (78) 

2.  Single  Wake  Case 

For  the  case  of  a  single  wake,  the  wake  weighting  function  reduces  to 

^  (79) 

where 

d,  =  coslirm  (80) 

and 

=  - e'^'* smlTTin  .  (81) 
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Plots  of  the  single  wake  lift  deficienc>  fi  ction  are  shown  in  Figures  12  through  15. 

The  propulsive  force  coeffic  :it  an  be  calculated  by  replacing  F  and  G  in 
equations  (34),  (35),  and  (36)  with  f  tnc  G\  For  pure  plunge  (a^  =  0) 

,  (82) 


for  pure  pitch  (h^  =  0) 


(83) 


and  for  coupled  pitch  and  plunge 


^4 


•"“A 


(84) 


Plots  of  the  propulsive  force  coefficier.t  ;;s  a  function  of  fi’equency  ratio  for  k  =  0.1234 
are  shown  in  Figures  16  through  18.  It  is  interesting  to  note  that  for  0.24  <  m  <  0.7 
the  propulsive  force  is  greater  than  tiit  "ICatzmayr  effect"  with  the  maximum  occurring 
just  prior  to  m  =  0.5.  Comparing  Figure ;  9,  10,  and  11  with  Figures  16,  17,  and  18 
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Figure  12.  Single  wake  lift  deficiency  function  as  a  function  of  wake  spacing  for 
m  ■  0. 
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Figure  13.  Single  wake  lift  deficiency  function  as  a  function  of  wake  spacing  for 
m  -  0.25. 
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Figure  15.  Single  wake  lift  deficiency  function  as  a  function  of  wake  spacing  for 
m  “  0.75. 
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Figure  17.  Propulsive  force  coefficient  in  pure  pitch  as  a  function  of  wake  spacing 
with  a  single  wake. 
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Figure  18.  Propulsive  force  coefficin  t  ii  coupled  pitch-plunge  as  a  function  of  wake 
spacing  with  a  single  wake. 


respectively,  the  maximum  propulsl\c  force  (at  approximately  m  *=  0.5)  for  a  single 
wake  is  always  greater  than  the  infirr*  n  imber  of  wakes.  TTie  reason  for  this  increase 
is  that  the  absolute  value  of  C  and  C'  a  e  never  greater  than  1.0,  but  C*  as  m 
approaches  0.5  exceeds  1.0  for  certait;  -.vake  spacings. 

When  m  *  0.5,  this  correspon  ic  to  the  wakes  being  180°  out  of  phase  as 
shown  in  Figure  19,  and  the  vortex  fn  in  the  wake  of  the  reference  blade  in  the  current 
revolution  is  directly  above  an  oppositely  spinning  vortex  from  the  single  wake  of  the 
previous  blade  or  revolution.  This  efftcic  it  use  of  the  vortices  in  the  wake  causes  C* 
to  exceed  1.0.  As  the  number  of  wakes  i  icreases  the  vortices  from  lower  wakes  are 
not  always  180°  out  of  phase  with  the  le'crence  blade.  When  m  =  0.5,  each  wake  is 
180°  out  of  phase  with  the  wake  above  a  id  below  it.  Thus,  each  wake  is  in  phase 


38 


with  the  wakes  ±2  wake  spacings  from  it.  Since  the  finite  wake  weighting  function 
contains  the  tenn  e  '*’"",  the  value  of  each  term  in  the  wake  weighting  function  will 
alternate  in  sign.  However,  since  the  wake  spacing  is  greater  for  each  subsequent 
wake,  the  absolute  value  of  each  term  is  not  as  large  as  the  first  term,  and  thus,  the 
effect  of  the  in-phase  wakes  are  not  as  strong.  As  the  wake  spacing  increases  this 
effect  diminishes. 

3.  Multiple  Wake  Case 

When  the  number  of  wakes  increases,  the  value  of  frequency  ratio  at  which 
the  maximum  propulsive  force  coefficient  occurs  changes  as  shown  in  Figures  20 
through  22  for  h  -  2.0.  This  change  is  analogous  to  the  estimation  of  a  square  wave 
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by  a  finite  number  of  terms  in  a  Fourier  series,  where  the  propulsive  force  coefficient 
determined  using  Loewy's  wake  weighting  function  is  considered  the  square  wave.  As 
Figures  16  through  18  and  20  through  22  show,  an  increase  in  either  wake  spacing  (h) 
the  number  of  wakes  {N)  will  decrease  the  maximum  propulsive  force  that  can  be 
obtained.  The  limit  of  this  decrease  will  be  Loewy's  case  for  a  hover. 
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Figure  22.  Propulsive  force  coefficient  in  coupled  pitch-plunge  as  a  function  of  the 
number  of  wakes  at  h  *=  2.0. 
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IV.  APPLICATIONS 


A.  EXTENSIONS  AND  COMPARISON  OF  THEORY 
1.  Forward  Airspeed 

Extension  of  the  finite  wake  theory  can  easily  be  made  to  include  forward 
airspeed.  Shipman  and  Wood  [Ref.  4]  expand  on  Loewy's  work  incorporating  forward 
airspeed  via  the  advance  ratio  (p)  into  the  rotary  wing  unsteady  aerodynamic  problem. 
The  definition  of  reduced  fi'equency  becomes 


k  = 


mb 

(1  +p)ft/ 


(85) 


The  Biot-Savart  Theorem  is  used  with  wakes  extending  out  to  infinity.  Since 
freestream  velocity  at  any  point  on  the  rotor  blade  is  a  function  of  both  blade  section 
radius  and  azimuth  position,  the  shed  vorticity  will  experience  a  build-up  and  decay  as 
azimuth  varies  from  0'  to  SbO”.  Shipman  and  Wood  developed  a  decay  function  to 
account  for  the  variation  in  vorticity  and  incorporate  this  function  and  the  advance 
ratio  into  the  integral  downwash  equation.  The  pressure  distribution  is  foimd  by 
applying  Sohngen's  inversion  formula  [Ref.  14]  to  the  downwash  equation.  The 
pressure  distribution  is  shown  to  be  in  the  same  form  as  Theodorsen  and  Loewy,  and  a 
modified  lift  deficiency  function  is  defined  to  account  for  the  differences.  Since  this 
modified  lift  deficiency  function  includes  Loewy's  wake  weighting  function  as  one  of 
its  terms,  the  finite  wake  weighting  function  could  be  directly  substituted  into  Shipman 
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and  Wood's  lift  deficiency  function.  If  the  segment  of  the  wake  weighting  function 
due  to  the  vorticity  decay  (AW)  is  modified  to  account  for  a  finite  number  of  wakes 
instead  of  an  infinite  number,  a  new  finite  wake  forward  airspeed  lift  deficiency 
function  can  be  defined,  and  the  remainder  of  the  problem  is  the  same  as  Shipman  and 
Wood. 

2.  Compressibility  Effects 

Adding  the  effects  of  compressibility  to  the  rotary  wing  unsteady  aerodynamic 
problem  was  first  accomplished  by  Jones  and  Rao  [Ref.  5].  Hammond  [Ref  6] 
provided  an  alternate  approach  to  the  problem  shortly  after  Jones  and  Rao,  but  only 
the  method  by  Jones  and  Rao  will  be  presented  here. 

Jones  and  Rao  expand  on  Loewy's  work  by  taking  the  known  solution  to  the 
fixed  wing  two-dimensional  unsteady  airfoil  theory  and  modifying  it  to  include  the 
effect  of  an  infinite  number  of  wakes  below  the  reference  airfoil  (or  rotor  disk). 

Euler's  equation  is  used  with  a  coordinate  change  to  determine  the  downwash  equation. 
An  interesting  result  is  found,  and  that  is  for  a  given  circulation,  the  downwash  due  to 
an  infinite  system  of  wakes  is  the  same  for  compressible  and  incompressible  flow.  In 
other  words,  the  circulation  generated  by  layers  of  shed  vorticity  beneath  an  airfoil  is 
independent  of  Mach  number.  This  fact  will  be  true  whether  there  is  an  infinite 
number  of  wakes  or  a  finite  number  of  wakes.  'Thus,  the  finite  wake  theory  can  be 
applied  here,  and  the  remainder  of  the  problem  is  the  same  as  Jones  and  Rao. 


43 


3.  NFS  Unsteady  Panel  Code 


The  NPS  Unsteady  Panel  Cr  i  v  as  written  by  Pang  [Ref.  7]  to  solve  the 
potential  flow  for  two  airfoils  execut ,  ^  u  isteady  motion  in  an  inviscid, 
incompressible  flow.  In  this  code,  tl  tw  •  airfoil  surfaces  are  approximated  by  a  large 
number  of  surface  elements,  and  a  ui.  or  i  source  distribution  and  vorticity  are  placed 
on  each  element.  The  source  streng  1  va’  es  from  element  to  element  ,  while  the 
vortex  strength  is  the  same  for  all  e.e  en  ;.  The  singularity  strengths  are  determined 
from  the  flow  tangency  condition  or  '  ot)  airfoil  sairfaces  and  the  two  Kutta 
conditions. 

Extensive  comparisons  of  th.-  co  le  with  Theodorsen  was  accomplished  by 
Riester  [Ref.  17],  and  the  code  (using  an  NACA  (XX)7  airfoil)  shows  very  good 
agreement  with  flat  plate  theory.  Con  pa  ison  of  this  code  with  Loewy  could  not  be 
accomplished  easily  since  it  would  rec  jire  a  very  large  number  of  airfoils  to  generate 
the  infinite  number  of  wakes.  Using  tl  •  finite  wake  theory,  the  code  could  easily  be 
compared  to  the  case  of  a  single  wake  A  comparison  for  pure  plunge  [Ref.  8]  with  h 
*  2.0,  k  *  0.1234,  and  hg  =  0.14  are  presented  in  Figure  23.  The  results  show  that  the 
NPS  Panel  Code  has  very  good  agree;  iiei  t  with  the  finite  wake  theory. 

B.  HIGHER  HARMONIC  CONTROI 

Higher  Harmonic  Control  (HHC)  is  in  active  vibration  control  system  for 
helicopters.  The  concept  of  HHC  is  to  a  ter  the  aerodynamic  loads  on  the  rotor  blades 
such  that  the  blade  response  is  reduced,  vhich  in  turn  reduces  the  vibratory  forces  and 
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Figure  23.  Comparison  of  propulsive  force  coefficient  in  pure  plunge  {h  -  2.0,  k  = 
0.1234,  and  -  0.14). 

moments  acting  at  the  hub  that  cause  vibration.  HHC  flight  test  results  [Ref.  15], 
under  a  joint  NASA-Army  sponsored  program  using  a  modified  OH-6A,  show  a 
reduction  of  airframe  vibrations  up  to  90  percent.  Basically,  HHC  is  an  electronic, 
computer-controlled  active  vibration  suppression  system  which  senses  and  cancels 
vibrations  in  a  helicopter  by  Q  per  revolution  feathering  or  pitch  motion  of  the  rotor 
blades,  where  Q  is  the  number  of  rotor  blades.*  On  the  OH-6A,  higher  harmonic 
blade  pitch  control  was  achieved  by  superimposing  4/rev  (or  4P)  swashplate  motion  on 
basic  cyclic  and  collective  control  inputs.  TTie  aircraft  was  flown  from  hover  to  100 
knots,  but  the  remainder  of  the  discussion  will  focus  on  the  hover  performance. 

*Noniially  for  helicopters,  N  is  used  to  represent  the  number  of  blades,  but  to 
remain  consistent  with  Oiapter  III  Q  will  be  used. 
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In  addition  to  reducing  vibration  levels  HHC  on  the  OH-6A  showed  the  potential 
for  decreased  helicopter  power  requirements.  Figures  24  and  25  show  hover  power 
required  for  three  different  schedules  of  HHC  open  loop  excitation.  What  is  of 
interest  here  is  that  the  amount  of  power  reduction  achieved  (on  the  order  of  10%  or 
20  horsepower)  appears  to  be  relatively  independent  of  the  type  of  excitation  applied. 
This  power  reduction  is  important,  especially  in  hover  where  power  required  is  large, 
because  of  the  benefit  of  lowering  fuel  consumption,  which  can  be  translated  into  a 
larger  payload  or  increased  mission  time. 

The  physics  behind  HHC  can  be  explained  by  unsteady  aerodynamics.  The  input 
by  HHC  to  the  rotor  blade  is  a  pitch  oscillation,  ±V3  degree  on  the  OH-6A.  From 
Figure  10  (infinite  wakes)  it  can  be  seen  that  the  propulsive  force  coefficient  is 
negative  (drag)  for  the  wake  spacing  of  the  OH-6A  (h  =  2.0)  regardless  of  the  phasing 
of  the  wakes  (m).  Thus,  it  might  be  concluded  that  pitch  oscillations  alone  cannot 
account  for  the  substantial  reduction  of  power  recorded  by  the  flight  tests.  If  it  is 
assumed  that  viscous  interaction  of  vortices  reduce  the  problem  to  a  finite  wake 
solution  as  shown  in  Figures  17  and  21,  the  propulsive  force  obtained  is  still  not 
sufficient  to  account  for  the  power  reduction  even  if  the  wakes  were  optimally  phased. 

As  stated  earlier,  HHC  provided  a  4P  pitch  input  in  both  collective  and  cyclic  were 
applied  at  the  blade  root.  The  collective  input  results  in  a  4P  oscillation  along  ihe 
blade,  while  the  cyclic  input  results  in  3P  and  5P  oscillations  along  the  blade.  From 
the  Southwell  plot  of  the  OH-6  A  [Ref.  16]  shown  in  Figure  26,  it  can  be  seen  that  the 
second  flapwise  and  third  flapwise  natural  frequencies  occur  near  3P  and  5P 


46 


Figure  24.  Effect  of  HHC  on  main  rotor  shaft  power  in  hover  for  selected  set  of  three 
open  loop  inputs. 

respectively.  During  the  course  of  a  revolution  of  a  blade  about  the  hub,  the  blade 
will  naturally  experience  flapping  motion  due  to  changes  in  angle  of  attack,  or  pitch 
changes,  as  the  blade  proceeds  around.  Since  3P  and  5P  flapping  frequencies  are 


Figure  25,  Effect  of  HHC  on  main  rotor  shaft  power  in  hover  for  selected  set  of  two 
open  loop  inputs. 


near  natural  frequencies,  small  3P  and  5P  pitch  inputs  will  result  in  much  greater 
flapping  outputs.  Thus,  the  unsteady  aerodynamics  of  HHC  for  the  OH-6A  is  coupled 
pitch-plunge  oscillations.  Figures  11,  18,  and  22  show  that  coupled  pitch-plunge 
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Figure  26.  Southwell  plot  of  the  OH -6 A  rotor  blade. 


oscillations  always  result  in  a  propul^i  e  >rce  regardless  of  the  number  of  wakes 
considered  or  phasing  of  those  wake  lat  is  of  greater  importance  is  that  if  the 
wakes  are  optimally  phased  together,  le  ]  ropulsive  force  is  an  order  of  magnitude 
higher  than  pure  plunge  or  pure  pile)  Tl  t  bottom  line  for  the  OH-6A  is  that  the 
inherent  natural  frequencies  of  the  bl  .long  with  4P  cyclic  inputs  result  in  coupled 
pitch-plunge  oscillations  which  yielc  .  pr  pulsive  force  large  enough  to  account  for  a 
substantial  power  reduction  for  the  1  c  co;  er. 
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V.  CONCLUSION 


The  unsteady  aerodynamics  of  rotary  wings  is  inherently  more  complicated 
because  of  the  layers  of  shed  vorticity,  or  wakes,  beneath  the  rotor  disk  that  interact 
with  the  rotor  blades.  The  Biot-Savart  Theorem  is  used  to  account  for  the  effects  of 
these  wakes,  and  a  modified  lift  deficiency  function  is  developed  for  use  with 
Theodorsen's  equations  of  motion  [Ref.  1].  Loewy  [Ref.  2]  develops  a  lift  deficiency 
function  for  the  case  of  an  infinite  number  of  wakes,  while  this  thesis  develops  a  lift 
deficiency  function  for  any  finite  number  of  wakes.  This  thesis  also  presents  the 
special  case  of  a  single  wake,  and  shows  that  for  small  wake  spacings  there  are  certain 
reduced  frequencies  where  the  lift  is  actually  enhanced.  In  other  words,  the  lift 
deficiency  function  becomes  a  lift  efficiency  function. 

The  equations  describing  the  propulsive  force  generated  by  an  oscillating  airfoil 
are  developed  by  Garrick  [Ref.  3].  These  equations  are  derived  from  Theodorsen's 
equations  of  motion,  and  included  the  lift  deficiency  function  as  a  parameter.  Since 
the  rotary  wing  lift  deficiency  functions  (Loewy's  and  the  finite  wake)  are  just 
modifications  of  Theodorsen's  work,  Garrick's  work  can  be  applied  to  rotary  wing 
urjjteady  aerodynamics.  When  either  Loewy's  lift  deficiency  function  or  the  finite 
wake  lift  deficiency  function  is  used,  the  propulsive  force  can  be  greatly  enhanced 
with  the  proper  phase  relationship  of  the  wakes.  The  use  of  coupled  pitch  and  plunge 
can  result  in  a  propulsive  force  which  is  an  order  of  magnitude  higher  than  either  pure 
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plunge  or  pure  pitch.  This  coupling  of  plunge  and  pitch  is  used  to  explain  the 
performance  benefits  reported  in  the  flight  test  data  for  the  OH-6A  equipped  with 
Higher  Harmonic  Control.  Coupled  pitch-plunge  oscillations  with  proper  phase 
relationships  of  the  wakes  yield  a  propulsive  force  large  enough  to  account  for  a 
substantial  power  reduction  for  the  helicopter. 
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APPENDIX  A 


I.  FLIGHT  OF  THE  HUMMINGBIRD 

For  thousands  of  years,  man  has  been  intrigued  by  the  flight  of  birds.  Early  Greek 
mythology  includes  a  tale  of  Daedalus  and  Icarus  who  made  wings  held  together  by 
wax,  and  escaped  their  imprisonment  on  the  island  of  Crete  by  flying  through  the  air 
like  birds.  Early  attempts  by  man  to  fly  were  usually  designed  aromid  imitations  of 
birds.  Most  of  these  early  attempts  met  disastrous  results  because  of  a  lack  of  clear 
understanding  of  the  aerodynamics  of  ornithological  flight,  and  that  is  tlie  wing  on  a 
bird  is  a  coupled  lift  and  propulsion  device.  In  fact,  the  Wright  brothers  first 
successful  flight  demonstrated  man's  limited  knowledge  of  aerodynamics  by 
decoupling  lift  and  propulsion  in  order  tc  achieve  flight.  To  this  day  lift  and 
propulsion  remain  decoupled. 

The  coupling  of  lift  and  propulsion  in  birds  can  be  explained  by  unsteady 
aerodynamics.  As  a  bird  flaps  its  wing  down,  the  trailing  edge  also  begins  to  rise 
(Ref.  18].  'Thus,  the  motion  of  the  wing  is  a  coupled  pitch-plunge  motion.  It  was 
shown  earlier  that  coupled  pitch-plunge  motion  yields  a  propulsive  force  that  is  an 
order  of  magnitude  higher  than  pure  plunge.  Since  most  birds  require  some  forward 
airspeed  before  being  capable  of  achieving  flight  and  the  rate  of  flapping  is  between  1 
to  5  beats  per  second,  the  analogy  to  unsteady  aerodynamics  can  be  described  by 
Garrick's  equation  for  coupled  pitch -plunge  using  Theodorsen'  lift  deficiency  function 
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(equation  (36)).  Birds  take  advantage  of  the  extra  propulsive  benefit  of  coupled  pitch- 
plunge  to  make  flight  appear  effortlessly. 

A  few  birds,  most  notably  the  hummingbirds,  are  capable  of  hovering.  In  hovering 
flight  the  hummingbird's  body  is  angled  upward  at  about  45®,  so  that  the  plane  of 
wing  flapping  is  approximately  horizontal.  As  the  wing  moves  up  (backward)  and 
down  (forward),  the  angle  of  attack,  or  p  tch,  changes  making  a  figure  eight  pattern 
that  is  on  its  side.  In  essence,  the  hummingbirds  wing  motion  is  similar  to  the  motion 
of  the  hands  of  a  swimmer  when  he  is  treading  water.  Once  again  the  hummingbird 
takes  advantage  of  a  coupled  pitch-plunge  motion  to  achieve  a  higher  propulsive  force. 

What  is  unique  about  the  hummingbird  is  the  rate  of  flapping,  or  the  oscillation 
frequency;  it  is  on  the  order  of  40  to  80  beats  per  second.  Welty  [Ref.  18]  states  that 
the  hovering  flight  of  the  hummingbird  is  grossly  inefficient  when  compared  to 
helicopter  rotor  blades,  but  this  is  overcome  by  an  increase  in  apparent  velocity 
created  by  trailing  currents  from  the  previous  reverse  stroke.  This  is  not  exactly 
correct.  'The  hummingbird's  wing  motion  is  similar  to  that  of  helicopter  rotor  blades. 
Both  are  cyclic  motions  of  changes  in  pitch  (angle  of  attack)  and  plunge  (flapping). 

'The  hummingbird  is  just  as  efficient  as  a  helicopter  -  both  have  a  large  power  required 
in  hover.  'The  apparent  increase  in  velocity  can  be  explained  by  trailing  vortices. 

Since  the  hummingbird's  rate  of  flapping  is  large  in  a  hover,  the  wake  created  by  a 
previous  stroke  will  not  dissipate  and  be  carried  downstream.  Each  stroke  of  the  wing 
always  has  layers  of  shed  vorticity  directly  below  it.  If  these  wakes  are  optimally 
phased  depending  on  the  number  of  wakes,  the  benefit  will  be  a  propulsive  force 
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greater  than  that  of  non-hovering  birds.  The  hummingbird  must  use  wakes  efficiently, 
or  it  will  not  hover.  Thus,  the  analogy  of  the  hummingbird  to  unsteady  aerodynamics 
can  be  described  by  Garrick's  equation  for  coupled  pitch-plunge  using  the  finite  wake 
lift  deficiency  function,  or  in  the  case  of  a  pure  hover  Loewy's  lift  deficiency  function. 
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